DNA Elasticity : Topology of Self- Avoidance 
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We present a theoretical treatment of DNA stretching and twisting experiments, in which we discuss 
global topological subtleties of self avoiding ribbons and provide an underlying justification for the 
worm like rod chain (WLRC) model proposed by Bouchiat and Mezard. Some theoretical points 
regarding the WLRC model are clarified: the "local writhe formula" and the use of an adjustable 
cutoff parameter to "regularise" the model. Our treatment brings out the precise relation between 
the worm like chain (WLC), the paraxial worm like chain (PWLC) and the WLRC models. We 
describe the phenomenon of "topological untwisting" and the resulting collapse of link sectors in 
the WLC model and note that this leads to a free energy profile periodic in the applied link. This 
periodicity disappears when one takes into account the topology of self avoidance or at large stretch 
forces (paraxial limit). We note that the difficult nonlocal notion of self avoidance can be replaced 
(in an approximation) by the simpler local notion of "south avoidance". This gives an explanation 
for the efficacy of the approach of Bouchiat and Mezard in explaining the "hat curves" using the 
WLRC model, which is a south avoiding model. We propose a new class of experiments to probe 
the continuous transition between the periodic and aperiodic behavior of the free energy. 



PACS numbers: 82.37.-j,36.20.-r,87.f 5.-V 
I. INTRODUCTION 

Motivated by experiments [1,2] in which single DNA 
molecules are stretched and twisted to measure their elas- 
tic properties, Bouchiat and Mezard (BM) [3] have pro- 
posed the Worm Like Rod Chain (WLRC) model, which 
gives a fair fit to the experimental data. However, several 
theoretical aspects of the WLRC model remain unclear, 
as evidenced by recent letters in PRL [4,5]. Points of 
dispute are BM's use of a local Fuller writhe formula (as 
opposed to the non-local Calugareanu- White formula [6] ) 
and the need for a phenomenologically introduced cutoff 
parameter, which has to be adjusted to fit the data. 

The theoretical issues raised by the experiments of 
[2] are surprisingly subtle. A reading of Bouchiat and 
Mezard's discussion of the "local writhe formula" [3] may 
give the misleading impression that the Calugareanu- 
White formula and the Fuller formula [7,8] are simply dif- 
ferent ways of expressing the same quantity, the writhe. 
In fact, there can be no local writhe formula. The writhe 
which appears in Calugareanu's theorem Lk = Tw + Wr 
is given by the non-local Calugareanu- White formula and 
jumps by two units when the curve is passed through 
itself. No local integral can replicate this behaviour. 
Bouchiat and Mezard use the "local writhe formula" un- 
der the assumption that the Euler angles are regular func- 
tions of the arc length parameter s. This assumption 
breaks down at the coordinate singularities of the Euler 
angles. So, Bouchiat and Mezard are effectively comput- 
ing the distribution of the "Fuller writhe" , which is not 
the same as the Calugareanu- White writhe. Neverthe- 
less, Bouchiat and Mezard obtain good agreement with 



experiment. The question remains: Why does it work? 
The purpose of this paper is to clarify these issues: we 
show that these questions are related to the topology of 
self avoidance. 

This paper is structured as follows: we first describe 
the experimental setup and summarize the experimen- 
tal data. We then explain the problem that we address, 
which concerns a global topological subtlety in the con- 
figurations of ribbons. In resolving the problem we point 
out the connection between the Worm Like Chain (WLC) 
[9], the Paraxial Worm Like Chain (PWLC) and the 
Worm Like Rod Chain (WLRC) models. We propose 
an experiment which brings out the relation between the 
models by continuously interpolating between them. We 
conclude the paper by comparing with previous work. 



II. EXPERIMENTS 

Micromanipulation techniques are now so sophisti- 
cated that experimenters can stretch and twist single 
DNA molecules to probe their elastic properties under 
a torsional constraint. In a typical experiment [2], one 
end of a single molecule of double stranded DNA is at- 
tached to a glass plate and the other to a magnetic bead. 
The glass plate is kept fixed and the bead is pulled by 
magnetic field gradients and rotated by magnetic fields. 
By such techniques the molecule is stretched and twisted 
by hundreds of turns and the extension of the molecule is 
measured as a function of the number of turns applied to 
the bead for a fixed force [2]. Another experiment mea- 
sures the torque-twist relation at fixed force [10] using a 
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slightly different experimental technique. The length of 
the DNA molecule is typically about 20 fj,m, its thickness 
2 nm and the bead is about 4.5/um in diameter. In prac- 
tice, this size of bead is adequate to prevent the molecule 
from untwisting by looping around the bead [4] . A typi- 
cal experimental plot is shown in Ref. [2]. These curves 
are easy to understand qualitatively: as one twists the 
molecule, its extension progressively decreases as one can 
guess by playing with a cord or ribbon. Further twist- 
ing leads to buckling and the formation of twisted braids 
or "plectonemes" which are familiar on telephone cords. 
Electron micrographs [11] of DNA show branched poly- 
meric structures which indicate the formation of plec- 
tonemes. 

To understand the experimental curves quantitatively 
more work is needed. Unlike telephone cords, DNA is 
seriously affected by thermal fluctuations and there are 
cntropic effects to be accounted for. Several papers have 
already treated this problem [3,4,12,13]. The high force 
regime, which is amenable to perturbation theory about 
the taut polymer configuration is well understood [12,13]. 
In the low force regime where self-avoidance effects are 
appreciable [14], a consensus is still lacking and this is 
the primary focus of this paper. 

III. THEORETICAL MODELS 

To define a theoretical model we need to specify the 
allowed configurations C of the polymer given the exper- 
imental constraints, write down a microscopic Hamilto- 
nian or energy E(C) for each configuration CeC and ex- 
press the partition function as a sum over configurations 
with Boltzmann weight. All experimentally accessible 
quantities can be got from the partition function. The 
problem therefore is to calculate the partition function 

Z(f 7 Lk) = J2e- E(C)/kBT (1) 

CeC 

where r is the vector separation between the ends of the 
molecule, Lk is the number of times the bead has been 
turned (which could be fractional) and the summation 
appearing in the expression for the partition function 
represents a sum over all allowed configurations of the 
polymer. E(C) is an energy functional which assigns an 
energy E(C) to each allowed configuration C. Specifica- 
tion of the allowed configurations and this energy func- 
tional defines the model completely. The next step is 
to "solve" the model, i.e deductively work out its experi- 
mental consequences. The objective is to confront theory 
with experiment and learn from the discrepancy as well 
as the agreement. Unfortunately, "solving" the model 
is not always a practical proposition even in idealised 
situations. Solvable theoretical models are therefore a 
valuable aid to understanding. 



It is often convenient to deal with conjugate vari- 
ables (B, f) instead of (Lk, r) and the Fourier-Laplace 
transform of the partition function. Z(B,f) = 
j dLkdrZ(Lk,r)e~ lBLk+ f ' r . We shall restrict ourselves 
to the limit of long polymers (i.e. for a polymer of con- 
tour length L and persistence length Lp, L/Lp — > oo). 
In this limit, the (B, /) and (Lk, r) ensembles are equiv- 
alent. This equivalence holds only in the limit of long 
polymers [15]. 

A twist storing polymer of length L is modelled [3,12] 
as a ribbon {x(s),u(s)} where s , < s < L is an ar- 
clcngth parameter along {x(s)}. u represents the "ribbon 
vector" and is required to satisfy u -t = where t = 4-x. 
The ribbon is described as a family of curves x(s)+eu(s), 
where x(s) is a curve and u(s) represents a slight (e) 
thickening of it along u(s) with e a small parameter. Let 
k(s) = \-^t\ be the curvature and £ls(s) = t ■ [u x j^u], 
the twist of the ribbon. In the models we discuss, the 
energy E(C) of a configuration C is given by 

E[C] = 1/2 / ds[AK 2 (s) + Cflj(s)], (2) 
Jo 

where A is the bending modulus and C the twist mod- 
ulus. Other terms can be added as in [12] with more 
parameters, but these are not necessary for our purposes. 

In order to compute the partition function Z(r,Lk), 
at fixed link Lk and extension r, we need to specify C, 
the allowed configurations of the polymer. This is what 
distinguishes the different models we now describe. In 
all the models we consider here, the configurations of the 
polymer are required to obey the following constraints: 
the polymer extends from x(0) = to x(L) = r, has 
fixed tangent vectors at the ends, t(0) — t(L) — z and 
the ribbon vectors at the ends are given by u(0) = x and 
u(L) — R(2nLk)x. R(2nLk) represents a rotation about 
the z axis through Lk turns. In addition, there may be 
further constraints which define the model and such con- 
straints can alter the topology of the configuration space 
and qualitatively alter the predictions of the theoretical 
model. For quantitative agreement with experimental 
data one needs to take into account the geometry and 
the statistical mechanics of the model. 

An important principle to bear in mind is that in Eq. 
(1), one should only sum over a single topological class 
[16] . Thermal agitation makes the polymer explore differ- 
ent configurations. But since these agitations only cause 
continuous changes, the polymer will remain in a sin- 
gle topological class. We should therefore sum over all 
configurations in a single topological class and not sum 
over distinct topological classes. Once this principle is 
understood and consistently applied, we find that the 
theoretical picture becomes much clearer. 

(a) Worm like chain (WLC): In this model no fur- 
ther constraints are imposed on x(s). Thus the poly- 
mer is allowed to intersect itself. While the set of 
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self-intersecting configurations may be of small measure 
in the configuration space C a , such configurations pro- 
foundly affect the topology of the configuration space. 
For, a WLC polymer can release link two units at a time 
by passing through itself [17,16,18]. This "topological 
untwisting" results in a collapse of link sectors. All link 
values differing by 2 are in the same topological sector 
and one has to sum Z(f, Lk) over Lk classes differing by 
2. This model has been treated in ( [16]) and leads to 
a partition function which is periodic in Lk with period 
two. This periodicity is clearly at variance with the ape- 
riodicity seen in the experimental curves [2] , where as we 
mentioned before, the extension progressively decreases 
as the bead is turned hundreds of times. The discrepancy 
between the WLC model and experimental data is clearly 
due to not taking into account the topological effects of 
self avoidance. 

However, there is a limit (the paraxial worm like chain 
(PWLC)) of the WLC model which yields an aperiodic 
free energy profile. In the high tension regime, (large 
force or \r\ comparable to the contour length L), the 
polymer is essentially straight between its ends and the 
tangent vector t only makes small deviations from the di- 
rection z of the applied force, which we call the north pole 
(of the sphere of tangent directions). In this regime, one 
can do perturbation theory and approximate the sphere 
of directions by the tangent plane at the north pole. This 
defines the 

(b) Paraxial worm like chain (PWLC): the tan- 
gent vector t = must be near z, the direction of the 
applied force i.e. t <~ z. 

The PWLC model has been treated in [12] and a simple 
analytic form for the writhe distribution is given in [13]. 
In this model, the polymer cannot release link by passing 
through itself since the high force prevents the molecule 
from looping back on itself to release link. However, this 
model only works in the limit of large forces (theoreti- 
cally, infinite) and does not address the low force regime, 
which is experimentally accessible. In order to prevent 
the polymer from releasing link at low forces by passing 
through itself, one would like to impose a condition pre- 
venting the polymer from intersecting itself. This defines 
the 

(c) Self- Avoiding Worm Like Chain (SAWLC): 

the configurations must be self avoiding: 

x(s) ? x(s'), (3) 

for s ^ s'. However, even this condition (3) is not suffi- 
cient to prevent topological untwisting. The experiments 
reported in [2] study the elastic properties of linear DNA 
molecules, not circular ones. So we need to address the 
issue of modelling open ribbons. Open ribbons again 
need careful handling because of topological untwisting: 
an open ribbon can untwist itself by two turns even in 
model (c) by going around its end and thus can change 



its link Lk by 2. Consequently, a configuration with a 
link Lk is in the same topological class as a configuration 
with link Lk + 2 and we once again get a collapse of Lk 
sectors leading to a periodic partition function [16]: 

Z{f,Lk + 2) = Z(f,Lk). 

In real experiments, there is an obstruction to releasing 
two units of link by going around the ends: the size of 
the magnetic bead is large enough to prevent topological 
untwisting of the polymer over the duration of the exper- 
iment. To model this obstruction to topological untwist- 
ing we close the open ribbon with a reference ribbon that 
goes through the bead, makes a large fixed circuit and re- 
turns to meet the polymer at the glass slide. We require 
that the reference ribbon and the real ribbon together 
form a non-self-intersecting closed ribbon. (Note that we 
can allow for a fractional link by permitting a discontinu- 
ity in the ribbon vector u at s = L: u(L — e) = Ru(L + e) 
where R represents a rotation about t and e is a small pa- 
rameter.) The condition of self avoidance keeps the real 
ribbon from going around its ends and effectively con- 
strains the applied link. It is understood, that when we 
sum over polymer configurations, we only sum over con- 
figurations in which the real polymer is changed. This is 
the "active part" of the polymer, as opposed to the "pas- 
sive" reference ribbon which is not summed over. The 
self avoiding worm like chain is defined by only allowing 
configurations which do not intersect themselves or the 
reference ribbon. This configuration space is called C c . 

As we mentioned before, this model is not tractable 
analytically. Ideally we would like to work with C c and 
sum over (the active parts of) all non self-intersecting 
closed ribbons extending from to r (and returning to 
along a reference curve) with fixed Lk. However, the 
constraint of self-avoidance is far from easy to handle! 
The constraint is non-local in the arc length parameter s 
and is hard to implement analytically [14]. This has led 
Bouchiat and Mezard [3] to propose a new model which 
they call the Worm like rod chain (WLRC). 

(d) Worm Like Rod Chain (WLRC) : is de- 
fined by imposing the additional constraint on the WLC 
configuration space C a that the tangent vector t must 
nowhere point south t ^ —z. Recall that the force is 
applied in the z direction, which we call north. This is 
also the direction in which the tangent vectors at the end 
are held fixed. Unlike the self avoidance constraint, the 
"south avoidance" constraint is local in the arc length 
parameter s. As a result, the model is amenable to semi- 
analytic treatment. The authors of [3] enforce the con- 
straint by using a repulsive potential at the south pole of 
the sphere of directions. The width of this potential is e, 
which serves as a cutoff. This model gives an aperiodic 
free energy profile and a very good fit to the experimen- 
tal data after the cutoff is suitably adjusted. While the 
WLRC model does agree with the data, its theoretical 
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significance is not clear to many workers in the field (see 
the recent exchanges in PRL [4,5] over the significance 
of the cutoff and the local writhe formula) . The purpose 
of this paper is to explain the nature of the relationship 
between the WLC, the WLRC and the PWLC models, 
explore different regimes of the elasticity of twist storing 
polymers and note that one can continuously interpolate 
between them in experimental situations. To summa- 
rize, the self avoiding Worm Like Chain (SAWLC (c)), 
the paraxial worm like chain (PWLC (b)) and the Worm 
Like Rod Chain (WLRC (d)) emerge by placing addi- 
tional constraints on the Worm Like Chain WLC model 
(a) [16]. One important qualitative feature that distin- 
guishes WLC from the other models is the fact that WLC 
predicts a periodic free energy whereas the other models 
predict aperiodic free energies. In the WLC model the 
system is delocalized over Lk sectors. This can be visual- 
ized in terms of a periodic multiple well potential where 
the thermal energy is high enough for a particle to visit 
all wells. From this point of view SAWLC, WLRC and 
PWLC are topologically and qualitatively similar and 
markedly different from the WLC model. 

The topological differences between these various mod- 
els can also be expressed in terms of the variable B 
conjugate to the applied link Lk variable. This conju- 
gate variable B, can be interpreted as a magnetic field 
[3,13,16]. In the a) WLC model B has the interpreta- 
tion of the field of a magnetic monopole located at the 
center of the sphere of tangent directions to the polymer 
[16] and the magnetic flux is quantized and this leads 
to a periodicity of the free energy. In the b) PWLC 
model, the excursions of the polymer are confined to 
the tangent plane at the north pole, and the magnetic 
field B is on a plane and therefore not quantised. As 
a result the free energy is not periodic in the link Lk. 




FIG. 1. The free energy of the WLC model has a period- 
icity (due to a collapse of Lk sectors) which can be removed 
by imposing self- avoidance or south-avoidance or large stretch 
forces. All these are different ways of separating the Lk sec- 
tors and preventing the "collapse of link sectors" . 



In the WLRC model, because of the constraint of south 
avoidance, B is defined on a sphere punctured at the 
south pole. Since the topology of the punctured sphere 
is identical (by stereographic projection) to that of the 
plane, there is no quantization condition just as in the 
PWLC model and again, this leads to an aperiodic free 
energy profile (see Fig.l). 

IV. WREATHE AND WRITHE 

In fact the right model to describe a real polymer un- 
der a torsional constraint would be c) the self-avoiding 
WLC (SAWLC) model which puts a constraint on the 
polymer passing through itself or the reference ribbon. 
However, as discussed earlier, such a model is hard to 
solve because of the non- locality of the constraint. We 
show below that the WLRC model captures some of the 
essential topological features of the SAWLC model and 
thus enables us to deal with a simpler and solvable model 
where the constraint reduces to a local one instead of a 
non-local one. Below, we show that i) the WLRC model 
captures the right qualitative behaviour of the Partition 
function. The WLRC partition function is aperiodic in 
the applied link unlike the WLC partition function [16] 
which is periodic, ii) Quantitatively, the WLRC partition 
function is a better approximation to the SAWLC par- 
tition function than the PWLC partition function and 
works even for low forces. 

We begin with the qualitative features, which are of a 
topological nature. The central quantity of interest [Eq.l] 
is the partition function Z{Lk) which is the sum over rib- 
bon configurations with fixed link [Lk). Calugarcanu's 
theorem [6] tells us that Lk = Tw + Wr, where the twist 
Tw = J Q 3 (s)ds is the integral of the local twist ^(s) 
along the curve. Writhe (Wr) is a quantity that only 
depends on the curve {x(s)} and not on the ribbon vec- 
tor {u(s)} [14]. The problem thus neatly splits into two 
parts. The link distribution is the convolution of the 
twist distribution and the writhe distribution [3,13]. The 
twist distribution is a quadratic path integral and is eas- 
ily evaluated. The problem that remains is to compute 
the writhe distribution. Thus the problem reduces from 
one defined on the space of ribbon configurations (C) to 
one defined on the space of ribbon backbones (C) or the 
space of curves. We can formally set the twist elastic 
constant to infinity (C = oo, the molecule is impossible 
to twist) so that Lk — Wr and then the Lk distribution 
is identical to the Wr distribution. 

The writhe is a non-local quantity defined only on 
closed non-self intersecting curves: Let s range over 
the entire length Lq of the closed ribbon (real ribbon 
+ reference ribbon) and let us consider x(s) to be a 
periodic function of s with period L . Let R(s,a) = 
x(s + a) — x(s). R(s,a) is non-vanishing for a ^ 0, L 
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and the unit vector R(s, a) is well-defined. It is easily 
checked that as a — > {0, L }, R — ► {t, — £} respectively. 
The Calugareanu- White writhe is given by [18-20] 



^ /"-£<o rL — 

Wcw =-J ds L da[ 



dR{s,a) x rfi?(s,ff) 1 _ ^ ^ 



(is 



G?<7 



The non-locality of Wciv makes it difficult to handle 
analytically. However, the key point to note is that vari- 
ations in Wcw are local [7]. Let C(r) be a family of 
non-self intersecting closed curves with writhe Wcw{t). 
Taking the r derivative of Eq. (4) we find that the re- 
sulting terms can be rearranged to give [21] 



dWcw 1 r° , t dt dt. ~ 

— ; = 7T~ <P ds\ — x — • t 

dr 2n J ds dr 



(5) 



which clearly has the interpretation of the rate at which 
t(s, t) sweeps out a solid angle in the space of directions. 
Note that Eq.(5) is a single integral [21] and therefore 
a local, additive quantity. Changes of writhe are local 
integrals and get contributions only from the active parts 
of the curve. Eq. 5 can be rewritten as: 



dWcw 1 rfO(r) 



dr 



2ir dr 



(6) 



where f2 is the solid angle enclosed by the oriented curve 
{i(s)\0 < s < L } on the unit sphere of tangent direc- 
tions as s goes from to L [14,16]. Note that £1 is only 
defined modulo 4ir: for a solid angle to the left of the 
oriented curve t(s) is equivalent to a solid angle (in — fi) 
to the right. dQ/dr is however well defined and local. 
Integrating Eq. (6) we arrive at [7]: 



Ztt 



2n 



(7) 



where n is an arbitrary integer. 

There is a quantity one can construct from the writhe 
which is well defined on all curves (not just simple ones) 

w(C) = exp[iirW C w(C)} = -exp[iSl/2] 

is a complex number of modulus unity which we call the 
wreathe. From the geometric phase point of view (see 
the analogy developed in [16]) the wreathe is a very nat- 
ural object to consider: it is simply the geometric phase 
of a spin half system in a cyclically varying magnetic 
field. When a curve is passed through itself, Wcw jumps 
by 2, but the wreathe is unchanged. We can therefore 
smoothly extend w to all closed curves including non- 
simple curves. The wreathe can be used to define a lo- 
cal quantity, the "Fuller writhe" [7] for curves which are 
nowhere south pointing. 

Let us define the "wreathe angular velocity" A T = 
—iw^ 1 ^- = ^dfl/dr as a "vector potential" on the space 
of curves. It is easily seen that A T is curl free. However, 



A T is not a gradient in the space of all curves: there exists 
no function W defined on all curves such that At is given 
by ^f. This follows because there exist closed circuits 
C(t) in the space of curves for which the integral § drA T 
is non zero. Such circuits enclose a nonzero topological 
flux and are link changing closed circuits (LCCCs). They 
describe the process of topological untwisting. 

Let us choose a fiducial curve C*, which goes from 
to r in a straight line, whose tangent vector is identically 
north pointing. Observe that all south avoiding curves 
are deformable to the fiducial curve C* via south avoiding 
curves. One simply deforms the tangent vector t(s) along 
the unique shorter geodesic connecting t(s) to the north 
pole. We now define the Fuller writhe as 



drA T - 1 



Writing the unit tangent vector as 
i = (sin cos 0, sin sin 0, cos#), J Q dr^ can be written 
as J ds^(l — cos 9) for all curves for which the tangent 
vector never points towards the south pole of the sphere 
of tangent directions. We can therefore write a "local 
writhe" on such curves which we call "south avoiding 
curves" : 



W F = ^- [ ds(l - cos (9)^ - 1. 
2tt J ds 



(8) 



While (8) is expressed in local co-ordinates on the sphere, 
it has a clear geometric meaning: 2tt(1 + Wf) is equal to 
the solid angle swept out by the unique shorter geodesic 
connecting the tangent vector t to the north pole [22,23]. 
This definition is explicitly not rotationally invariant, 
since it uses a fixed fiducial curve C* and singles out a 
preferred direction. 

This definition of "Fuller writhe" is motivated by a 
theorem of Fuller [7], which gives the CW writhe of a 
curve which is deformable to the fiducial curve by a fam- 
ily of curves which is self avoiding and south avoiding. 
For such curves, the Fuller writhe agrees with the CW 
writhe. However, in general, the "Fuller writhe" is not 
equal to the CW writhe, which is what appears in the 
Calugareanu- White formula. In fact, the definition of 
Fuller writhe extends easily to all curves which are south 
avoiding, since Eq. (8) does not require any more condi- 
tions. 

We remarked earlier that A T is not a gradient on the 
space of all curves. Only in certain restricted classes 
of curves can it be expressed as a gradient. For exam- 
ple, A is a gradient in the space of self avoiding curves 



. dWcw 
(It 



. It is also a gradient in the space of south 



avoiding curves A T = n^ff- [24]. 

To summarise, Wcw is defined on all self-avoiding 
curves, Wf on all south avoiding curves. On curves 
which arc deformable to the fiducial curve C* through 
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south and self avoiding curves the two notions agree [7] 
(Wcw = Wf). When a curve passes through itself, 
Wcw jumps by two and when the tangent vector to 
a curve swings through the south pole, Wf jumps by 
two units. The writhe is a real number which has both 
geometric and topological information. The topological 
part is the integer part of and the geometric part 
is the fractional part of The geometric part is com- 
pletely captured by wreathe but the topology is lost since 
wreathe is insensitive to changes in writhe by 2 units. 
From the definitions it is clear that on curves where these 
quantities are well defined, 

.vj- 1 dw dWcw dW F /n . 

—i = = (9) 

ir dr dr dr 

so changes in writhe are the same whether measured by 
y^cw, VVf or w. Note that w is well-defined on all curves. 

Let us now apply these general topological ideas to un- 
derstand the configuration space of the polymer ribbon. 
As stated earlier, we need to sum over a single topolog- 
ical sector of the configuration space. For a ribbon, this 
means summing over a single knot class K of the ribbon 
backbone C and a single link class Lk of the ribbon. So, 
in fact the true partition function ZK(r,Lk) would de- 
pend not only on the extension and the link, but also the 
knot class of the ribbon backbone. Needless to say, this 
is a hopelessly intractable problem. 

The constraint of self avoidance is a non-local con- 
straint in s and very hard to handle analytically. Even in 
the pure bend model [25,26] one cannot handle this con- 
straint and so one just gives up self avoidance and sums 
over all configurations. This is of course an approxima- 
tion, but it seems to work well. Instead of evaluating (1) 
with C equal to simple curves in a single knot class K 
(for example, the trivial knot class) : 

z * = E e ~ mc) ( 10 ) 

CeK 

we are effectively summing over knot classes and com- 
puting Z = Y] K Zx since we are unable to impose self- 
avoidance. The resulting analysis still gives a reasonable 
account of the data [1,25]. This suggests that the con- 
tribution from the non-trivial knot classes may not be 
significant. In fact, in the presence of a force /, there 
is an energy cost fLp leading to a suppression e~t Lp of 
the probability of forming knots [27] . 

In the case of twisting polymers, giving up self- 
avoidance has a more serious consequence. Giving up 
self avoidance leads not only to an identification of knot 
sectors, but also link sectors separated by two units. This 
is due to LCCCs. These are closed circuits in the space 
of curves, which when lifted up to the space of ribbons 
by continuity become open and lead to an identification 
of link sectors. We show below that LCCCs must pass 
through both south pointing and self intersecting curves. 



More precisely, our main result is: in any closed circuit, 
the number of signed self crossings is equal to the number 
of signed south crossings. 

Proof: The proof uses the fact that the wreathe is 
defined on all curves. Consider a closed circuit C(r) 
(0 < t < 1) in C a which starts from C* and returns 
to it: C(0) = C(l) = C„. If we now compute the wreathe 
(which is defined on all curves) we find that as r varies 
from to 1, < r < 1, the wreathe w(t) describes a mo- 
tion on the unit circle and returns to its starting point. 
The number of times w(t) winds around the unit circle 
is given by 

k= — [ drw^ 1 dw I dr = — <f> drA T 

7T Jo * J 

Using Eq. (9) we find that k measures i) the number 
of times (counted with sign) the polymer passes through 
itself (i.e. k = 1/2 ^ dr^ff^-). ii) the number of times 
(counted with sign) the polymer tangent vector swings 
through the south pole (i.e. k = 1/2 ^ dr^f-). 

The above discussion can be "lifted" in a sense (made 
precise in the appendix) to the space of ribbons which 
is the true configuration space. Starting with a closed 
circuit C(t) and an initial ribbon C» whose base curve is 
C„. One can lift the circuit C(r) in the space of curves, 
by continuity to the space of ribbons. However, closed 
circuits in the space of curves may be open in the space 
of ribbons! This is the well known anholonomy effect. 
Our main result can then be re-expressed on the space of 
ribbons: two ribbons based on the same ribbon backbone 
C* are homotopic as self avoiding ribbons if and only if 
they are homotopic as south avoiding ribbons. 

Proof: The argument above also shows that 2k is 
equal to the net change in the Lk class of a ribbon if 
its backbone is continuously deformed along C(r) (since 
Lk = Wr in our analysis). 

If we specialise to self avoiding curves C c then k = 
AWcw — 0- If w e work with south avoiding curves Cd, 
then k = AWf = 0. This proves that self and south 
avoidance present the same topological obstruction to re- 
lease of link Lk. This provides a formal justification for 
Bouchiat and Mezard's work [3] where they impose a 
south avoidance constraint in place of a self avoidance 
constraint. 

k also has the interpretation of the "quantized mag- 
netic flux" passing through the loop {C(t);0 < t < 1}. 
We can write k as 

r 1 r L ° di di If 
k = l/(2n)<b dr & ds[— x — ] • t = - & drA T (11) 
Jo Jo ds dr 7T J 

which measures the "topological flux" passing through C. 
k is also the degree of the map t(r, s) from the torus {0 < 
s<£o>0<T<l}to the sphere of tangent directions. 
Clearly if the flux is non zero, then the degree is non 
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zero and the tangent vector must point south somewhere. 
Removing the south pole [28] (as Bouchiat and Mezard [3] 
do) forces the degree to be zero and prevents link collapse. 

Our arguments above show that the topology of self 
avoidance is captured by south avoidance. This is a qual- 
itative feature related to topology. At this qualitative 
level, the excised point on the sphere of tangent direc- 
tions need not be the south pole (recall that the force 
direction is taken to be the north pole), but can be any 
direction h on the sphere. However, in order to achieve 
quantitative success in explaining the experimental data, 
we show below that it is advantageous to choose h along 
the south pole. Ideally we would like to compute the par- 
tition function of the SAWLC model to compare with ex- 
periment. The replacement of SAWLC by WLRC seems 
to give a reasonable account of the data. Why does it 
work? To see why, consider the region 1Z in the space of 
curves which can be reached from the fiducial curve C* 
without ever self intersecting or pointing south. Let us 
notice the following: i) The region 1Z has a finite measure 
in the space of curves, ii) All over 1Z, Wcw = Wf. iii) 
The region 1Z dominates the partition function. More 
precisely, regions outside 1Z are suppressed by a factor 
e -f eff L P _ Here j eff = f-B 2 /(4L P ) is an effective force 
[12,13] that takes into account the competition between 
the stretch and the twist. Consequently, the WLRC par- 
tition function is a good approximation to the SAWLC 
partition function over a range of parameter space. 

Let us successively consider the regimes of low, in- 
termediate and high twist. Quantitatively, this is mea- 
sured by Jeff- (i) For f e ffLp large, the tangent vec- 
tor makes small excursions about the north pole and 
< 9 2 >= f e ffLp is small. In this regime, one can 
approximate the sphere by its tangent plane at the north 
pole and as explained earlier, this is the regime of the 
PWLC model, (ii) For < 9 2 > of order l,(/ e //L P of 
order 1), this is not a good approximation, since the ge- 
ometry of the sphere is not well approximated by the pla- 
nar geometry. However, in this regime, the WLRC works 
well since the suppression factor e~^" ffLp is small. Since 
the WLRC retains the geometry of the sphere (changing 
only the topology, by removing a point), the WLRC is an 
improvement on the PWLC model, (iii) At very low (or 
negative) f e ff, in the SAWLC model, the polymer ac- 
commodates writhe by winding around itself. A similar 
behaviour occurs in the WLRC model where the tangent 
vector is mostly near the south pole. In both models a 
cutoff set by the thickness of the polymer (see below) is 
needed to regularise the model. 

Had we chosen to exclude some other direction instead 
of the south pole, the suppression factor would have been 
e feffL P cose ^ wnere q i s p i ar an gle of the excluded 

point. Clearly, the advantageous choice is 9 = tt i.e., 
the south pole of the sphere. The choice 9 = n leads 
to an axially symmetric quantum mechanical problem, 
for which the Hamiltonian can be solved. This is exactly 



what Bouchiat and Mezard choose in their WLRC model. 

Need for a cutoff: In summing over all south avoid- 
ing paths one finds that the problem of computing Eq. 
(1) is ill defined. Recall that Sd, the set of south point- 
ing paths have been removed from C a , the WLC con- 
figuration space to get the WLRC configuration space 
Cd = C a — Sd- The problem is that near Sd there are 
paths (points in Cd) with arbitrarily large writhe and 
vanishingly small energy. This problem is due to paths 
which wind around the south pole in tiny circles accu- 
mulating large writhe at zero energy cost [29]. 

This problem does not occur at large forces. At large 
forces there is a large energy cost preventing the tangent 
vector from visiting the south pole. Indeed, in the parax- 
ial limit there is no pathology [4,12,13]. As one lowers the 
force, the energy cost for accommodating writhe tends to 
zero near Sd and these dominate the sum in Eq. (1). In 
order to get sensible results from the WLRC model one 
has to impose a cutoff: the tangent vector is excluded 
from a small circle 9 = 9 C = n — e around the south pole. 
In more detail, we see that for paths winding around the 
south pole the energy per unit length E/L ~ e 2 while 
the writhe per unit length Wr/L ~ 2 — e 2 . This implies 
that in the limit e — > any amount of writhe can be 
accommodated at zero energy cost. Maggs and Rossetto 
[4] claim that this pathology of the WLRC model is an 
artefact of the local Fuller writhe formulation. However, 
we notice that the same pathology afflicts the SAWLC 
model (which involves a non local notion of writhe, the 
CW writhe) as well. This can be seen as follows. Con- 
sider a polymer twined on a cylindrical rod [30] like a 
garden climber twines around a pole. As explained in 
Fuller [7], for a plectoneme of pitch angle a, (in the limit 
that a goes to ir/2 the helix is almost a straight line) the 
curvature goes as k = (d/2)~ 1 cos 2 a, while the writhe 
per unit length Wcw /I goes as sina, where d is the di- 
ameter of the rod. (See Fuller, figure 1 and its caption 
in the 1971 paper). From the geometry cosa = nd/(l), 
where I is the length of the polymer. When a goes to 
7r/2, the writhe (per unit length) goes to a constant and 
the energy per unit length goes as d 2 . Note that we are 
holding the length of the polymer constant as we take 
the limit a going to tt/2. 

We thus come to the conclusion that both SAWLC (in- 
volving the non local CW writhe) and WLRC (involving 
the local Fuller writhe) suffer from the same pathology 
which can be cured by introducing a cutoff parameter 
which has the same origin ( the finite thickness of the 
polymer) in both formulations. In fact Bouchiat and 
Mezard's WLRC model is a simple analytically tractable 
model which captures the essential qualitative physics of 
the SAWLC model. 
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V. TUNNELLING 

So far we have used topological ideas to split up the 
configuration space into disconnected sectors. However, 
the difference between topology and energetics can be 
blurred by activated processes [4]. To appreciate this 
point, consider a Brownian particle in a double well po- 
tential. If the barrier height is not too large, the par- 
ticle will randomly visit both wells during the observa- 
tion time. If one increases the barrier height, the acti- 
vated process is exponentially suppressed. In the limit 
that the barrier height goes to infinity, we may say that 
the two wells belong to distinct topological sectors and 
sum over only one of them. The transition from one 
regime to the other is measured by the Kramers' time 
TKramcrs ~ exp(V/£;T), where V is the height of the bar- 
rier and kT the temperature. If T Kr amors >> ^observation, 
then we would sum over a single well, but if the reverse 
is true, then we would sum over both wells. By tun- 
ing V, the barrier between the wells, we can continu- 
ously interpolate between the single well and the dou- 
ble well descriptions. The picture can be extended to 
include multiple wells, for example, a potential profile 
V(x) = —VqCosx which has multiple minima at x = 2mr. 
In the polymer context, the different Lk sectors corre- 
spond to different minima of this potential. Tunnelling 
between different minima is an activated process and the 
classical solution that dominates this process is an "in- 
stanton" [31]. The physics involved here is very similar 
to the tunnelling between different topological sectors of 
vacuum classical configurations in QCD. The "0 vacuaa" 
that emerge in QCD are delocalised over all the topolog- 
ical sectors, just as the WLC polymer is delocalised over 
link sectors. Tunnelling between different minima is ruled 
out in the models (b), (c) and (d) but permitted in the 
model (a). In order to interpolate between these mod- 
els theoretically we can soften the cutoff in the WLRC 
model by putting a potential of height Vq and width e. As 
Vo — > oo we recover the WLRC model and as Vo — > we 
recover the WLC model. Vq serves as a natural param- 
eter in interpolating continuously between these models. 
In the next section we show how one can experimentally 
explore the continuous transition between these models 
using the catalytic effect of enzymes. We will see how 
addition of Topoisomerase II enzymes effectively lowers 
the energy barrier Vq permitting the molecule to pass 
through itself and thus enabling the polymer to be delo- 
calised over distinct link sectors. 

VI. PROPOSED EXPERIMENTS 

We now explore the experimental realization of the 
WLC model (model (a)). The process of "topological 
untwisting" can be made possible by using the enzyme 



Topoisomerase II to permit the DNA molecule to pass 
through itself [32,33]. (Similar experiments are referred 
to in [34].) As the concentration of enzyme c is increased 
we would expect a transition from aperiodic behavior 
to periodic WLC behavior. A low concentration would 
correspond to a large energy barrier to passing through 
itself and a high concentration lowers the energy bar- 
rier. We may expect the energy barrier Vq to go as 
c _1 . A similar effect is also expected to happen as a 
function of force /. At large forces (or extensions), the 
energy barrier to looping back and passing through it- 
self is prohibitive. As one lowers the force, one expects 
again a transition from the aperiodic behavior noticed in 
a paraxial worm like chain (PWLC) model Rcf.[ [12], 
[13]) to the periodic WLC behavior [16]. The height 
of the barrier can be worked out in the WLRC model 
( model (d)). It is given by Vq + 2/ where Vo is the 
repulsive potential at the south pole of the sphere of 
tangent directions and / is the applied stretching force. 
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FIG. 2. The theoretically expected "phase diagram" in the 
force-concentration (/ — c) plane, c is the concentration of 
the enzyme Topoisomerase II, which permits self crossing. 

Plotting the lines of constant Kramer's time gives a qual- 
itative phase diagram in the force-concentration (/ — c) 
plane (see Fig. 2). At high forces and low concentrations 
the free energy is aperiodic. At low forces and high con- 
centrations the free energy is periodic in the link. 

During the process of replication the DNA molecule 
undergoes supercoiling and it needs to unwind to release 
its stress. The viscous cell environment offers resistance 
to this process. The unwinding of these supercoiled struc- 
tures takes place via Topoisomerase enzymes which cut 
the polymer and helps it to release its stress. This real 
biological context is where the WLC periodic free energy 
finds a natural application. These effects can be studied 
under controlled circumstances in single molecule exper- 
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iments. 

VII. CONCLUSION 

In this paper we have provided a formal justification 
for the WLRC model in terms of the topological effects 
of self-avoidance. Our main result is that: In any closed 
circuit of curves, the number of signed self-crossings is 
equal to the number of signed south-crossings. Or cquiv- 
alently, in the language of ribbons, two ribbons based on 
the same ribbon backbone C* are homotopic as self avoid- 
ing ribbons if and only if they are homotopic as south 
avoiding ribbons. 

Replacing the non-local notion of self-avoidance by the 
local notion of south avoidance results in the analyti- 
cally tractable WLRC model. The Calugareanu- White 
formula for the writhe is explicitly nonlocal (the formula 
is expressed as a double integral) and no local formula 
valid on all simple curves exists. In Fuller's treatment 
[7] a theorem is proved that the writhe difference for two 
closed non self-intersecting space curves is given by a lo- 
cal formula under certain conditions [7,8]. Taking the 
special case that the reference curve has constant tangent 
vector z, we arrive at the formula [Eq. (8)] for non self- 
intersecting (i.e simple) curves which are nowhere south 
pointing. Note, however, that the restriction of non self- 
intersection can be relaxed since it nowhere appears in 
the Fuller writhe formula. Our treatment which uses the 
wreathe as the starting point has a larger domain of va- 
lidity since wreathe is well defined on all curves. Our 
main result provides the theoretical framework and jus- 
tification for the calculational scheme employed by [3] in 
the WLRC model. 

We have focussed on long polymers to keep this dis- 
cussion simple and because the experimental situation 
permits it. Needless to say, the general topological dis- 
cussion applies to short polymers as well. In particular 
consider a very rigid polymer (L p — > oo, L/L p — > 0) with 
the tangent vector clamped at both ends to z. This sys- 
tem can be modelled as a paraxial worm like chain (even 
if the applied force is zero) and according to the general 
discussion, will have an aperiodic free energy. In this 
case, the energy barrier to looping around is given by the 
elastic energy of the stiff polymer and goes as A/L. 

We have emphasized the need for introducing a regu- 
larizing cutoff in a theoretical analysis of this problem. It 
has been suggested [4] that the divergence in the writhe 
distribution seen in Ref. [3] is an artefact of the Fuller 
formulation of the writhe and will not be present in a 
CW formulation. We have shown that the CW writhe 
also suffers from the same pathology. One sees [7] that it 
costs (almost) no energy to accommodate any amount of 
writhe. As a result, large values of writhe are possible. 
This pathology is only cured by considering the physi- 
cal thickness of DNA (about 2 nm), which results in an 



energy cost for writhe and suppresses large writhe fluctu- 
ations. Writhe fluctuations can be suppressed either by 
increasing the tension on the molecule (PWLC model) 
or by taking into account the physical thickness of the 
molecule (WLRC model). 

Ribbons homotopic as self avoiding ribbons are also 
homotopic as south avoiding ribbons. Note that the con- 
verse is not true: knots cannot be undone if one imposes 
self avoidance. But one can remain within the set Cd 
of south avoiding curves and undo any knot. Consider 
projecting the knot onto the x — y plane and by small 
movements of the tangent vector to the polymer undo the 
knot by passing the ribbon through itself since there is 
no self avoidance constraint. When we sum over all south 
avoiding ribbons we are automatically summing over all 
knot classes of the central curve. This is an approxima- 
tion valid only when the contribution of the nontrivial 
knot sectors is negligible. 

In future it would be interesting to address the theo- 
retical issues that stem from considering distinct knot 
classes. The recent class of experiments [35] probing 
the dynamics of complex knots on single molecules of 
DNA would perhaps enable us to connect the theoret- 
ically challenging energetic and topological aspects of 
knots [27,17] to experiments. One can also experimen- 
tally explore the low force stretching regime to explore 
deviation from the pure bend WLC model of Ref. [25] to 
probe the effect of knotting of a DNA molecule. 

As an offshoot of the observation of the topological dis- 
tinction between various models of twisting polymers we 
are led to a new class of experimentally testable predic- 
tions. We hope this work will generate interest amongst 
experimentalists to test these predictions. 
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VIII. APPENDIX: RIBBONS, CURVES AND 
FIBRES 

To clearly put across the topological point at issue here 
we use some mathematical notation. Let C be the config- 
uration space of all closed ribbons and C the space of all 
curves. There is a natural map P : C — > C, which maps 
each ribbon C to its backbone curve C. The ribbons based 
on C arc classified upto homotopy by their link class Lk. 
The structure (C, C, P) constitutes a fibre bundle. Given 
a closed circuit in C, which starts from and returns to C* 
and a starting ribbon C* , one can continuously "lift" the 
circuit to the space of ribbons. However, a closed circuit 
in C may lift to an open circuit in C. The initial and 
final ribbons have the same backbone curve, but belong 
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to different link classes. We refer to such circuits as link 
changing closed circuits (LCCCs) in the text. 

Let C* be a fiducial ribbon whose backbone is C* . Our 
sum in Eq. (1) extends only over the component con- 
nected to C*. We call this component 7Tq(C). In order 
to perform the sum correctly we need to understand the 
structure of 7r (C) and 7Tq(C) correctly in the different 
models. 

We start with the model (c) the SAWLC model, which 
is the one of experimental interest. 7To(C c ) is labelled by 
K the knot class of the ribbon backbone C and Lk, the 
link class of the ribbon based on C. (Strictly speaking 
we should write [K] and [Lk] , the equivalence classes be- 
ing denoted by the square brackets, but we do not do so 
here). The Link is an integer and counts the Gauss link- 
ing number of the two edges of the closed ribbon. Thus 
7i"o(C c ) = (K,Lk). 7Tq(C c ), the component connected to 
C* consists of all ribbons in (if*, Lfc»), those in the knot 
class of C* (the unknot class) and which have the same 
link as C„. 

(a) WLC model: If one gives up self avoidance and 
works with C a , then the structure of 7To(C a ) is quite dif- 
ferent from C c . There are no knot classes since any knot 
can be undone by self crossings. The even link classes get 
identified with one another and similarly the odd link 
classes. This is due to the presence of Link Changing 
Closed Circuits (LCCCs) in C a . These curves are closed 
circuits in C a . If one lifts them up by continuity to C a , 
the space of ribbons, we find that the continuous lift of 
a closed curve in C a may be open in C a . The LCCCs 
can be recognised by the fact that they enclose a nonzero 
topological flux k = 1/n § drA T ■ The final ribbon has a 
link which differs from the initial one by 2k. Thus, the 
link classes of like parity are identified with each other 
and 7To(C a ) has just two components. 7Tg(C a ) includes 
all the link sectors which differ from C* by 2k. This is 
what we refer to as the "collapse of link sectors" . The 
configuration space C a has a vastly different topological 
structure from C c , the space of interest. This is where 
the south avoiding WRLC model (d) comes in. Knots 
can be undone in Cd but (as we see below) links cannot. 
So 7r (Cd) = Z, where Z is the set of integers, the link 
Lk. The claim that south avoidance captures the topo- 
logical effects of self avoidance is based on the following 
observation. Consider a LCCC with flux k. This provides 
a continuous deformation in C a , between ribbons (based 
on the same backbone), which differ in Lk by 2k. It is 
shown in the text that these LCCCs pass through self in- 
tersecting curves and south pointing curves k times (both 
counted with sign). Imposing self avoidance eliminates 
these LCCS. Alternatively we can permit self intersec- 
tions and impose south avoidance. This also has the same 
effect of eliminating LCCCs and preventing the collapse 
of link sectors, which occurs in C a . This proves that self 
and south avoidance both present the same topological 
obstruction to link release. Mathematically, the link class 



bundle over C a is nontrivial. One can locally trivialise 
this bundle by omitting points (in fact sets of measure 
zero) from C a . Restricting the bundle to Cd and C c 
gives two different trivialisations of the same bundle. 
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